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ABSTRACT 

Rapid variability on a time scale much faster than the light-crossing time of the central 
supermassive black hole has been seen in TeV emission from the blazar PKS 2155-304. The 
most plausible explanation of this puzzling observation is that the radiating fluid in the rel- 
ativistic jet is divided into a large number of sub-regions which move in random directions 
with relative Lorentz factors w 7'. The random motions introduce new relativistic effects, 
over and above those due to the overall mean bulk Lorentz factor Ft of the jet. We consider 
two versions of this "jets in a jet" model. In the first, the "subjets" model, stationary regions 
in the mean jet frame emit relativistic subjets that produce the observed radiation. The vari- 
ability time scale is determined by the size of the sub-regions in the mean jet frame. This 
model, which is motivated by magnetic reconnection, has great difficulty explaining the ob- 
servations in PKS 2155-304. In the alternate "turbulence" model, various sub-regions move 
relativistically in random directions and the variability time scale is determined by the size of 
these regions in their own comoving frames. This model fits the data much more comfortably. 
Details such as what generates the turbulent motion, how particles are heated, and what the 
radiation process is, remain to be worked out. We consider collisions between TeV photons 
emitted from different sub-regions and find that, in both the subjets and turbulence models, 
the mean bulk Lorentz factor Ft, of the jet needs to be > 25 to avoid the pair catastrophe. 

Key words: black hole physics - magnetic reconnection - relativistic processes - turbulence 
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1 INTRODUCTION 

On July 28, 2006, the high-frequency peaked BL Lac source 
PKS 2155-304 had a strong flare in the TeV band, with an av- 
erage flux that was mo re than 10 times large r than the typical 
flux seen at other times l lAharonian et al.ll2007h . During this flare, 
which lasted for about an hour, rapid variability on a time scale 
« 300 s was observed. With a galaxy bulge luminosity Mr = 
—24A JKotilainen et al.ll998h . the expected mass of the black hole 
(BH) in the nucleus of the galaxy is Mbh ~ 1 - 2 x 10^ Mq 
JBettoni et al.l l2003h . This corresponds to a gravitational radius 
Rg = 2GM/c^ « 3 — 6 X 10^'* cm, and a corresponding time 
scale Tg — Rg/c~l — 2x 10'* s. Amazingly, Tg is nearly two 
orders of magnitude longer than the variability time. 

It is generally accepted that TeV emission in blazars comes 
from relativistic jets tha t are pointed towards the observer 
JHinton & HofmannI 120091) . Opacity arguments suggest that the 
bulk Lore ntz factor of the jet in PKS 2155-304 must be very large. 
Lb > 50 teegelman et al.ll2008l) . At first sight, it might appear that 
this large Lorentz factor will also explain the rapid variability ob- 
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served in the source. However, while relativity can certainly cause 
the observed variability time to be shorter than R/c, where R is 
the size of the emitting region, there is no simple way for the vari- 
ability time to be shorter than the gravitational time scale Tg of the 
central engine. In PKS 2155-304, the variability is faster by a factor 
of several tens. 



The only way to understand the rapid variability in PKS 2155- 
304 is if one of the following conditions holds: (i) the entire source 
(engine and emitting region) moves towards the observer with a 
Lorentz factor ~ 50, or (ii) the emitting region alone moves rapidly 
towards the observer, and the variability is caused by some local in- 
stability in the radiating gas which is insensitive to any time scale 
associated with the central BH engine, or (iii) the supermassive BH 
is ~ 50 times less massive than we estimate. The first option is ob- 
viously impossible. We do not expect a 10^ Mq BH to move with a 
Lorentz factor of 50. The third option is also unlikely since the BH 
mass estimate is obt ained via the well-know n (BH mass)— (bulge 
luminosity) relation JMagorrian et al.ll 19981) , which has an uncer- 
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tainty of no more than a factor of a fevo We are thus led to the sec- 
ond option, but even this is highly non-trivial. The size of the emit- 
ting region is surely at least as large as Rg, which means, given the 
short variability time scale, that only a small fraction of the emitting 
volume can be involved in the variability. Yet, the large amplitude 
of the observed fluctuations suggests that the whole emitting region 
must contribute to the flare. There is an inherent paradox in these 
conflicting indications. 

Building on ideas developed earlier to explain variability 
in gamma-ray bursts (GRBs, see iLvutikov & Blandfora l2003l : 



lLvutikov|[200a ; lNaravan & Kumaill2bo9l:lLazaret alj|2009l). Gian- 
nios et al. (2009; see also lOiannios et alJuOlOl : iNalewaiko et al.l 
120111) proposed a "jets in a jet" model as a possible explanation 
of the TeV flare observed in PKS 2155-304. According to this 
model, a relativistic jet moves with an overall bulk Lorentz fac- 
tor Fj, towards the observer, but the TeV emission is produced 
in emitting regions that themseves move relativistically with re- 
spect t o the mean frame of the jet. As shown in the context of 
GRBs JNaravan & Kumaj|2009l : iLazar et al.ll2009h . such a model 
naturally produces large amplitude variability on time scales much 
shorter than the light-crossing time of the emitting volume. 

In this work we explore the "jets in a jet" model and obtain 
several constraints that must be satisfi ed for this model t o oper - 
ate. Our analys i s foll ows the work of iNaravan & Kumag ( 120091) 
and lLazar et al.l 1 120091) on the equivalent problem for GRBs. In §2 
we summarize the relevant observations of PKS 2155-304. In §3 
we discuss the overall issue of time scales in relativistic jets. We 
then describe the "jets in a jet" model and its two basic variants 
— "subjets" (§3.1), "turbulence" (§3.2) — and analyze the condi- 
tions under which each can fit the observations. Both variants of the 
model involve a large number of emitting regions moving in differ- 
ent directions. Photons emitted in one region can interact with those 
emitted in another, and if the optical depth is sufficiently large, col- 
liding high energy photons will pair produce and make the source 
opaque. We examine the optical depth to this pr ocess in §5 and 
consid er how the constraint Ft > 50 obtained bv lBegelman et al.l 
( 120080 is modified in the "jets in a jet" scenario. In §6 we summa- 
rize our results and consider the implications. Additional details are 
given in two Appendices. 



2 OBSERVATIONS 

The TeV flare under consideration in PKS 2155-304 was al- 
ready in progress wh en observations began on July 28, 2006 
dAharonian et al.ll2007l) . It lasted for about Tobs ^ 4000 s, after 
which the flux fell to a much lower value. For the analysis in this 
paper, we need to know the total duration T of the flare. In Ap- 
pendix A, we use a Bayesian analysis to estimate the probability 
distribution of T. We estimate the median duration of the flare to be 
Tinedian ~ 2Tobs ~ 8000 s, and the mean duration to be Tmcan = 
ln(rn,ax/TobB) Tobs ^ (3 - 10)Tob. « 12000 - 40000 s. For 
the present analysis, we choose T « 20000 s (conesponding to 
Tmax ~ a week). 

The second quan tity we need is the variability time scale. 
lAharonian et alj ( 120071) found very rapid rise times r^ ~ 100 — 
200 s and longer decay times Td ~ 200 — 600 s (see their Table 
1). They also defined a doubling time T2, and state that the fastest 



^ An exception would arise, of course, if the sy stem consists of a bin ary 
black hole and the smaller BH produces the flare jVolpe & Riegenl201 If) . 



T2 = 224 ± 60 s and the average T2 = 330 ± 40 s. Based on 
this information, we choose the characte ristic variability time to be 
St ^ 300 s (as in lBegelman et alj|2008h . The ratio of the total du- 
ration of the flare to the variability time is an important quantitjQ 
For the above choices, we find 



r 

Tt 



70, 



T ^ 20000 s, 



<5tRi300s. 



(1) 



During the 4000 s of observations of the flare in PKS 2155- 
304, there were 5 major pulses in the TeV emission. Scaling this 
number to the total duration T, we estimate the total number of 
pulses in the flare to be A^p « 25. From this we estimate the duty 
cycle ^ of the pulsed emission: 

^ = ^ ^ 0.4, N, « 25. (2) 

The numerical values given in equations l[T]l and (|2} are for our 
fiducial estimates of T and St. As a very extreme case, we also 
sometimes consider T = Tobs ~ 4000 s, St ~ 600 s, which give 
T/St Ri 7, C ~ 0.8. 

As mentioned earlier, the mass of the BH in PKS 2155-304 is 
estimated to be Msn ~ 1 — 2 x 10^ Mq, so its gravitational time 
scale is 

Mbh 



Tg ^ lOOOOAfg s, 



Mq 



109 Mo 



1 



(3) 



In the following we use a canonical value of Afg = 2, such that 
Tg matches our estimate of the flare duration T. However, there 
is no direct measurement of the BH mass, so we also occasionally 
consider a lower value: Mg = 0.5. 



3 THE MODEL 

Naively, one might think that relativistic motion can cause the ob- 
served variability time of a source to be arbitrarily small and so 
rapid variability can always be explained. However, the situation is 
more complicated. 

Consider a smooth homogeneous jet moving towards the ob- 
server with a Lorentz factor Ft. Let the jet be at a distance R from 
the central engine and have a radial width AR as measured in the 
"lab" frame, i.e., the frame of the host galaxy. The emission is 
beamed into an angle l/Ft, so a distant observer receives radia- 
tion from only a region of lateral size ~ R/Tt. The angular time 
scale tang, the time delay between the center and edge of the visi- 
ble patch, is ~ R/2crl. Since the jet is assumed to be smooth, any 
observed variability can be no faster than this|j 

In the comoving frame of the jet, the size of a causally con- 
nected region is R' ~ R/Tt, where here and throughout the paper 
we use a prime to distinguish length and time scales in the comov- 
ing frame from those in the lab frame. The radial width A_R' of the 
shell of material in the jet cannot be smaller than this size. (If it 
starts out smaller, it will quickly expand to this size as a result of 
internal pressure gradients.) Thus, two photons emitted simultane- 
ously in the fluid frame from the front and back of the shell will 



^ This q uantity is sometimes called the source variab ility V in the GRB 
literature iSari & Piranlll997l : lNaravan & Kumaj2009h 
■^ Throughout this paper, for simplicity, we ignore the cosmological time 
dilation factor (1 + 2:) = 1.116 between the host galaxy of PKS 2155- 
304 and the observer. This factor will be more important for a source at a 
higher redshift. In that case, time scales in the observer frame will need to 
be multiplied by (1 -|- z). 
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reach the observer with a time difference /!\R' /2cVb > R/2cr\. 
This radial time scale irad is of the same order as, or longer than, 
the angular time scale iang. 

Remarkably, R/2cVl is also the observed time difference be- 
tween two photons emitted by a single fluid element, one emitted 
when the fluid is at a distance R from the engine and the other 
emitted at a distance 2R. This dynamical time scale fdyn is rele- 
vant, for instance, if the jet slows down as a result of interacting 
with a surrounding medium. It is also relevant for a freely stream- 
ing jet, since it is the time scale on which the energetic particles in 
the jet are cooled by adiabatic expansion. 

For large values of Fb, all three time scales described above 
are much shorter than the Newtonian time scale R/c that one would 
compute based on the linear size of the source. Thus relativistic mo- 
tion is indeed capable of producing very short time scales compared 
to the dimensions of the emitting region. 

However, there is a fourth time scale. From a simple light- 
crossing argument applied to the power source of the jet, the ra- 
dial width of the emitting shell of material cannot be smaller than 
the size Rg of the central engine. Thus, the shell width in the lab 
frame must satisfy A_R ^ Rg, or equivalently, in the jet frame 
AR' ^ TiRg. This additional constraint means that the observed 
variability time scale cannot be shorter than Tg. However, in PKS 
2155-304, we have Tg ~ 20000 s, whereas the variability time is 
much shorter, St ~ 300 s. This is the problem we are faced with. 

In order to solve the above problem, we have to give up the 
assumption of a smooth jet and must allow the radiation observed 
at any given time to be emitted from a tiny region of the source with 
a size I' <^ R' . However, R' is the size of a causally connected re- 
gion. How can a significant amount of energy, as needed to produce 
the observed large amplitude pulses of TeV emission, be squeezed 
into a region much smaller than R'l The answer is (Narayan & 
Kumar 2009; Lazar et al. 2009; Giannios et al. 2009; see Appendix 
B for a Table comparing the notations used in these papers with 
that in the present work) to make the additional assumption that 
the radiating region under consideration moves relativistically with 
respect to the comoving frame of the jet. This causes the radiation 
from a tiny source region to be beamed into a narrow cone, thereby 
amplifying the observed luminosity without enhancing the energy 
requirement of the emittin g region. This is the ke y idea of the "jets 
in a jet" model. Note that lGhisellini et alj ( l2009f) suggest an alter- 
native way of reducing the size of the emitting region, which we do 
not consider in this paper. 

In models of magnetically accelerated jets, during the acceler- 
ation phase the radi us and Lorentz factor of th e jet generally scale 
as _R ~ ^bRg (e.g.. lTchekhovskov et alj2008l) . This gives a causal 
scale R' = R/Vt ~ TbRg- If acceleration ceases and the radiation 
is produced when the jet is in a coasting phase, then R! > V^Rg- 
To allow for both possibilities, we write: 



R — fc^bRg 



R' = ^ 

1 b 



fcViRg 



(4) 



where the numerical factor fc^i allows for coasting. In the fire- 
ball model of GRBs, during acceleration we have R ~ TbRg and 
it would appear that we could have /c < 1. However, energy dissi- 
pation in a baryon-loaded fireball typically occurs ony at a distance 



R 



or equivalently, idy 



icnginc, where t. 



engine 



is the characteristic time scale of the central engine dSari & PiranI 
1 19971) . Clearly fengine ^ Tg, and thus once again we obtain equa- 
tion lO with /c > 1. 

Events such as th e bright flare observed in PKS 2155-304 by 
lAharonian et alj 1 120071) are clearly rare. Most of the time the source 



is much fainter. To model this behavior we assume that the source 
has more or less steady low-level jet activity, but occasionally goes 
through short periods of time when the engine power becomes very 
much larger. We are concerned with the properties of these bursts 
of more energetic activity; one of these bursts presumably produced 
the bright flare seen in PKS 2155-304. Let icngino be the duration 
that the source spends in an energetic state. Clearly ^engine J? Tg. 



In addition, if ie 



is very small, radial spreading will cause the 



shell to expand to a width ~ -R/F^, as discussed earlier We thus 
write the radial width of the emitting region as 



AR — max ( ct, 



R 



R 



i-ongine, p2 I — /t* p2 — Jd, -^ — JdJcRg, 



R' 

i b 



AR' = U 



(5) 
where the engine duration factor fd ^ 1- The overall duration of 
the flaring activity is determined by AR, since the angular spread- 
ing time and dynamical time are always shorter than or equal to the 
radial time. Hence 



^^ - f HL - f f T 

C FfcC 



(6) 



In PKS 2155-304, T is roughly equal to Tg for our fiducial choice 
of the BH mass {Mg — 2). Thus the product fdfc must be of order 
unity. Since fd and fc are individually larger than or equal to unity, 
this implies that each is of order unity. If we assume a lower mass 
for the BH, e.g., AIc) = 0.5, then fdfc ~ 4 and we have some 
freedom in choosing the values of fd and fc . 

As explained above, the key idea of the "jets in a jet" model 
is that the radiating fluid in the jet is sub-divided into a number of 
independent volumes, each moving relativistically with respect to 
the mean frame of the jet. In order to explain the rapid variabil- 
ity, it is necessary for each sub-volume to be much smaller than 
the causality scale R' , or the engine-related scale AR' . The break- 
ing up of the jet into kinematically distinct sub-volumes cannot be 
related to the engine since the length scale involved is too small. 
Rather, it must result from a local instability of some sort. Such 
a situation can arise naturally in a highly magnetized outflow, for 
instance through magnetic reconnection or MHD turbulence. How- 
ever, in the following, we do not assume anything specific about 
the nature of the outflow or what causes the instability. 

There are tw o main variants of the "jets in a jet" model 
JLazar et al.ll2009f) . The distinction is whether the emission dura- 
tion of individual pulses in the flare is determined by physics in the 
comoving frame of the jet or in a frame moving relative to the jet. 

3.1 Subjets and reconnection 

In the model described by Giannios et al. (2009; see also Lyutikov 
2006, Lazar et al. 2009 in the context of GRBs), which we denote 
hereafter as the "subjets" model, magnetic field reconnection cells 
arise sporadically within the strongly magnetized jet fluid. Each re- 
connection event leads to the ejection of twin subjets of relativistic 
plasma with a typical Lorentz factor 7^ as measured in the mean 
frame of the jet. The subjets emit the observed TeV emission. A 
single pulse in the observed TeV lightcurve corresponds to a single 
subjet. Hence the observed duration of a pulse is equal to the time 
taken for the completion of a reconnection event as measured in the 
jet frame, divided by F;, (to transform to the observer frame). While 
this model is strongly motivated by magnetic reconnection, there in 
nothing in the following discussion that depends on the specific de- 
tails of this mechanism. Hence the model is valid for any process 
in which a local intability produces relativistic subjets. 

Each reconnection event dissipates the magnetic energy in a 
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certain characteristic volume in the jet frame. Let I' be the typical 
length scale of this volume, and let /3' = v' /c be the typical speed 
with which the magnetic energy flows into the central reconnection 
zone. The duration of the reconnection event, as measured in the jet 
frame, is 5t' ~ I' / (3' c, and so the observed variability time is: 



St: 



VbP'c 



(7) 



Clearly, if a subjet changes its direction by more than 1/7^ before 
consuming all the energy within the reconnection volume then the 
observed pulse will be shorteQ In this case, we simply redefine /' 
such that the subjet direction is constant to within I/7J during the 
time 5t' = I' / P' c. Making use of equation ((SJ, the ratio of the flare 
duration to the variability time is: 



T 
St 



.M'^. 



(8) 



For a given observer, the region from which radiation can be 
seen has a size ~ R' perpendicular to the line-of-sight and a size 
~ AR' — fdR' (measured in the jet frame) along the line-of-sight. 
Since each independent reconnection region has a volume ~ I''', 
the total number of subjets within the observed volume is 



ntot 



2/d 



2a,3 



flP 



(9) 



The factor of 2 is because each reconnection site produces two sub- 
jets moving in opposite directions. Purely from relativistic beam- 
ing, each subjet would illuminate a solid angle ~ '^ Hf in the jet 
frame. Since there might be additional beam broadening due to in- 
trinsic velocity fluctuations within the subjet, we write the solid 
angle illuminated by one subjet as 0,'^ = i^fj/'yf with /_, ^ 1. 
Assuming that subjets from different reconnection regions are un- 
correlated and are oriented randomly, a given observer receives ra- 
diation from a fraction ~ Q!j /47r of the subjets. Each visible subjet 
produces one pulse in the observed lightcurve. We thus estimate the 
duty cycle in this model to be: 



^ 



/. 



T 



fj f R' 



27f ./!/?'' \StJ 2/3' y-y'l' 



(10) 



Using our fiducial numbers for PKS 2155-304, viz., T/5t 
70 and ^ « 0.4, equation JlOl l gives 



fj 



: 6000. 



(11) 



We showed earlier that fd ~ 1, and by definition we have fj ^ 1. 
In addition, current understanding of r elativistic reconnection sug- 
gests that /?' ~ 0.1 l lLvubarskvll200^ . Writing /3' = O.l^fl^i, we 
thus find 



7;^2500/;/Vd"'/3-f^'. 



(12) 



This i s an extremely large value, particularly when we recall that 
JLvuba rskv 20 0^ 7' ~ -y/q, where a is the magnetization param- 
eter JKennel & Coronitill984h . We require a to be truly enormous, 
which leads to other problems (see below). Even if we set /?' = 1 
ij3'_i — 10), which is unlikely, we obtain 7^ ~ 80. The subjets 
thus need to move highly relativistically with respect to one other 
and with respect to the mean jet frame. Also, equation ^ indicates 

^ Correspondingly, the shape of the pulse in the light curve will no longer 
be determined by the onset and decline of the reconnection event, but will 
be determined by the motion of the jet. 



that there must be ntot ~ 7 x 10 l3'Si independent subjets, which 
is again very extreme. 

Extreme values of the parameters, e.g., T ~ 4000 s (which 
requires a lower BH mass Alg < 0.4 to maintain T ^ Tg), 
St ~ 600 s (the longest time scale consistent with the observa- 
tions), improve the situation somewhat. With /3' = 0.1, this gives 
7j « 170, Titot ~ 7 X 10^, still rather extreme, while with /3' = 1, 
we obtain reasonable values: 7^ ~ 5, ntot ~ 700. Alternatively, 
we could assume that the BH has a smaller mass than our fiducial 
value Mg = 2. For instance, if we take Mg — 0.5 and assume 
fc = 1, we obtain fd — 4. If we further set fj = 1, /?' — 1, then 



7, 



20 and ntot ~ 4 x 10*. 



An additional issue is that, in order for relativistic reconnec- 
tion to be energetically efficient, we need the inflowing magnetic 
field lines on the two sides of the reconnection reg ion to be aligned 
to within an angle ~ l/2y/a JLvubarskvl 120051) . Recalling that 
7j ~ \/o" and that the model requires a large 7^, this poses a serious 
problem. 



3.2 Relativistic turbulence 

Another scenario, which was advocated by Narayan & Kumar 
(2009; discussed further by Lazar et al. 2009) and which we de- 
note the "turbulence" model, is one in which relativistic turbulence 
is generated, possible because of some MHD instability in the jet 
fluid. As a result, blobs of fluid move in random directions with a 
typical Lorentz factor 7^ as measured in the mean jet frame. Each 
blob is roughly spherical in its own frame, which means it has a 
transverse size I' and longitudinal size /'/7i 1" 'h^ J^t frame. At 
any instant, the radiation from a given blob is focused into a solid 
angle ~ t^ H'l in the jet frame. However, over time the beam ori- 
entation might wander, s o we write the solid angle illuminated by 
a blob as flj = -k f, /'7'^. lNaravan & Kumari ( 120090 suggested that 
the velocity vector of a blob might wander by about a radian during 
the time a blob radiates, which corresponds to /j ~ 7^. However, 
more general situa tions, including /j « 1, can also be considered 
JLazar et al.ll2009l) . 

The duration of a pulse as measured in the jet framqf] is given 
by the longitudinal size of a blob, i.e., St' ~ I' j^'jC. Transforming 
to the observer frame, the observed pulse duration is: 



St: 



I' 



Fb'y'c 



(13) 



The ratio of the duration of the flaring event to the duration of a 
single pulse is then: 



T ,R!_ 

St ~ ^"'^^ V ■ 



(14) 



The total number of emitting regions within the observed vol- 
ume is (there is no additional factor of 2 here since each blob radi- 
ates into a single beam): 



, 'R'Y 1 fT 

ntot ~ Id {-fT I = 



I' J fh? \st) ' 



and the duty cycle is: 



fj (TV ^ fo f R' 



47f/j v^*; ' 47^ w/ 



(15) 



(16) 



Compared to the subjets model, we see that /3' — >■ 7', with an 



^ See lLazar et al.l I2009h for a discussion of additional relevant time scales. 
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additional factor of 2 because of the different number of beams per 
blob in the two models. 

Generally, we find that the turbulence model has a much eas- 
ier time satisfying the observational constraints. For example, with 
T/5t Ri 70, ^ ^ 0.4, fj = 1, fd = 1, we obtain -y'j ^ 5 and 
Ttt ot ~ 3000, which are qui te reasonable. If fj ~ 7^ (as suggested 
bv lNaravan & Kumadl2009l) . then 7^ « 7 and Utot = 800, which 
is again acceptable. 

The huge difference in the predictions of the subjets and tur- 
bulence models can be traced to differences in how the variability 
time scale is related to the blob size I' in the two models. In the 
subjets model we have 5t' ~ I' / j3'c, i.e., it is the crossing time of 
a region of size I' , as measured in the mean jet frame, at speed /3'c. 
In the turbulence model, on the other hand, we have 5t' ~ I' j'^'jC, 
i.e., it is the crossing time of a region of size I' , as measured in the 
frame of the moving blob, at speed c. The factor /?' ~ 0.1 arises in 
the subjets model because we believe reconn ection is limited to a 
speed substantially below c JLvubarskvlbOOSr) . There is no equiva- 
lent factor in the turbulence model, mostly by assumption, since we 
do not have as complete a physical picture of this model as we do 
for the subjets model. The second factor 7^ is because of another 
key difference between the two models. In the subjets model, the 
time scale for a reconnection event is determined by physics in the 
mean frame of the jet, since we assume that the reconnection cell 
is at rest in this frame. In the turbulence model, on the other hand, 
everything is determined by physics in the comoving frame of the 
blob. The net factor of 7j//3' between the two models is quite a 
large number and this leads to drastic differences in their predic- 
tions. 



4 OPACITY LIMITS 



As discussed bv lCavallo & ReesI jl978f) . if a source of non-thermal 
high energy radiation is spatially more compact than a certain limit, 
then photon-photon collisions will be very frequent and there will 
be copious electron-positron pair production. The source will then 
be opaque to its own high energy radiation and will be unable to 
produce the obser ved non-thermal s pectru m. Assuming a spatially 
homogeneous jet, ISegelman et alj ( 12008) estimated that the pair 
production optical depth in PKS 2155-304 is 



2 X 10'"i46i-7nnr 



■461300-'- 6 1 



(17) 



where L46 is the (isotropic equivalent) luminosity of the jet in units 
of 10*® ergs^^ and tsoo is the variability time in units of 300 s[j 
Setting L46 ~ fsoo ~ 1, the requirement r-y^ < 1 gives V^ > 50. 
In the "jets in a jet" model, the radiation received by the ob- 
server during any given pulse in the light curve comes from a single 
subjet or blob in the source. Since the emitting region moves with 
respect to the observer with a net Lorentz factor ~ ^h'^'j, one might 
be tempted to replace Fj, by F67j in equation {Yl\ to thereby obtain 
ri,7j > 50. This would loosen the constraint on F(, by a large fac- 
tor (at least 5, and potentially several tens). However, the argument 
is incorrect[j 



^ The expo nenent on Fj, in eq. 1171 depends on the spectrum of the source. 
IPiranI 419991) gives 4 -|- 2q instead of 6, where a is the spectral index (Fd oc 
jy—a^ The exponent 6 thus corresponds to a = 1, which is a reasonable 
v alue for blazars and GR B s . 

■^ iGiannios et all )200g|) derived a bound Fj > 9 for PKS 2155-304 using 
what appears to be, apart from a number of details, essentially this incorrect 
argument. 



The result given in l ll7b corresponds to a single radiating re- 
gion moving relativistically towards the observer. It describes the 
optical depth for the radiation to escape from its own local emis- 
sion region. Since there are many emission regions in the "jets in 
a jet" model, even after a beam escapes from its original blob, it 
is likely to encounter other beams of radiation on its way towards 
the observer. In order for the radiation to reach the observer, we re- 
quire the net optical depth due to all these encounters to be small. 
To calculate the corresponding optical depth, it is most convenient 
to work in the mean jet frame (the primed frame in our notation). 

Let us focus on the volume of the jet from which the observer 
can receive radiation. This volume has a size _R' x i?' in the two 
transverse directions and A_R' = /d-R' in the radial direction. Each 
fluid element in this volume will radiate roughly for a dynamical 
time idyn ~ H! jc. Multiplying this time by c to convert it to a 
length, the total four-volume that is visible to the observer is 



K 



(4) 



/di?" 



(18) 



Consider first the subjets model. Each subjet illuminates a 
conical volume with a solid angle fji^/^f, and so the average 
three-volume of the cone is {fjT^/Z'^f){r^), where r is the length 
of the cone inside the reference volume. For r distributed uni- 
formly between and Ft! , the mean three-volume of a subjet is 
{fjTv/12'yj^)R''^. The time for which a subjet shines (as measured 
in the jet frame) is I' / j3' c. Multiplying by this factor, and also by the 
number of subjets ntot given in equation l|9j, the total four-volume 
occupied by all subjets is 



K 



(4) 



.fi^ 



subjctB ~ 67f /J/3'4 



^"'(li) 



(19) 



Dividing by V^^^^^ and making use of equation dlOt . the fractional 
four-volume occupied by subjets is 



(4) 



1/ , . -JT 

p subjets r^^ /- 

Jsubjets = ~7774i ^ "o"^' 



V, 



(4) 



(20) 



i.e., it is roughly equal to the duty cycle of the observed light curve, 
^ ~ 0.4. This result is not surprising. 

A similar calculation can be done for the turbulence model. As 
before, the solid angle of each beam is fji^/^f, but the duration of 
the beam in the jet frame is ~ I' j^'jC. Repeating the same steps as 
above (using eas. ll5|[T6t . we find 



K 



(4) 



/jTT 



127: 



/6 f 2 
Id 



R 



"''It. 



ftur 



yd) 
turbulence 



K 



(4) 



k. 



(21) 



(22) 



The result is the same as for the subjets model. 

Before proceeding, let us define a "reference jet model" which 
consists of a homogeneous jet with a comoving volume R' x R' x 
fdR', moving with a Lorentz factor Fj towards the observer. By 
construction, the total duration of the observed high energy radi- 
ation from this hypothetical jet is the same as in PKS 2155-304: 
T ~ fdR' /Yi,c (see eq.[6](. We assume that the jet produces the 
same mean luminosity as that observed in PKS 2155-304. However, 
being homogeneous, it cannot produce the observed rapid variabil- 
ity. 

When applying equation Ml\ to the reference jet model, 
we must set the variability time equal to R' /TbC — T/fd ~ 
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(20000//d) s. As a result, the optical depth and limiting Lorentz 
factor become (for L46 ~ 1) 



reference jet model 



:3xlo«/d^,-^ 



^6>25/y^ 

(23) 
The larger size of the emitting volu me, compared to the model con- 
sidered by ISegelman et alj 1 120081) where the size was constrained 
by the observed variability time of 300 s, leads to a factor of 2 re- 
duction in the minimum bulk Lorentz factor Fj,. 

Consider now the "jets in a jet" model. In both the subjets 
and the turbulence versions of this model, the radiation of the criss- 
crossing beams occupies a fraction ~ ^ of the total four-volume 
Kotai- Since the "jets in a jet" model produces the same average 
luminosity as the reference jet model, the mean number density of 
photons in the jet frame must be the same. The only difference is 
that, in the "jets in a jet" model, a fraction ^ of the volume is occu- 
pied by radiation, and in these regions the local number density is 
higher than average by a factor ^^^. Any given beam of radiation 
will intersect many other beams on its way ouij. A fraction ^ of 
its path will be through other beams, each of which will on aver- 
age have a photon number density a factor ^"^ larger than in the 
reference jet model, and the remainder of its path will be through 
radiation-free regions. The net result is that each escaping photon 
will on average interact with exactly the same number of high en- 
ergy photons as a photon in the reference jet model. Thus, the opti- 
cal depth to pair production is the same in both models, i.e.. 



"jets in a jet" model : 



SxioVdr^-', 



r, > 25/]/^ 

(24) 
Note that, while the optical depth is the same in both the ref- 
erence jet model and the "jets in a jet," the latter model has the 
distinction of being able to produce the rapidly varying light curve 
observed in PKS 2155-304. Interestingly, the model does this with- 
out suffering any penalty in the pair production opacity. In fact, the 
"jets in a jet" model has a less restrictive limit on the bulk Lorentz 
factor of th e jet, Tf, > 25, compare d to the original limit. Lb > 50, 
obtained bv lSegelman et alj ( l2008h . 



5 IMPLICATIONS AND CONCLUSIONS 

The very rapid variability o bserved in the TeV flux of PKS 2155- 
304 JAharonian et akluOOTD , which is faster by a factor of 50 than 
the gravitational time scale of the central BH, is a remarkable puz- 
zle. Normally, even assuming a relativistic outflow, one does not 
expect to see variability faster than the gravitational time of the 
driving engine. PKS 2155-304 violates this expectation by a large 
factor. A possible way out — apparently the only way out — is 
to invoke some version of the "jets in a jet" model. This model 
was discussed earlier (though this name was not used) in con- 
nection with the variability of G RBs (Lvutikov & Blandf oral2003l ; 
lLvutikovll200^ : iNaravan & Kum ar 2009; Lazar et al. 200§)). The 
sam e idea has been show n to explain the variability in PKS 2155- 
304 JGiannios et alj200^ . 

In the "jets in a jet" model the fluctuating pulses in the ob- 
served radiation are not produced by the entire jet but by small 
sub-regions within the jet. These sub-regions move relativistically 
with respect to the mean frame of the jet. As a result each sub- 
region produces a beam of radiation that is focused tightly in the 



° This would not be true if J <g 1, but we are working in a different limit 
where ^ ~ 1. 



direction of its motion. Only a few sub-regions radiate towards the 
observer, but these appear anomalously bright as a result of rela- 
tivistic beaming. Thus, the model explains both the large apparent 
luminosity and the rapid variability, without violating any energy 
or light-crossing-time constraints. 

We have considered in this paper two versions of the "jets in 
a jet" model. The first version, the "subjets" model, assumes that 
something like magnetic reconnection takes place in many differ- 
ent sub-regions within the jet fluid and that each reconnection site 
ejects twin relativistic subjets in opposite directions along the local 
recombining magnetic field. This model is characteristic of a class 
of models in which the time scale of pulses in the observed light 
curve is determined by a process that is at rest in the frame of the 
jet. The role of the relativistic subjet is merely to provide a luminos- 
ity boost through beaming. The second version, the "turbulence" 
model, assumes that some instability in the jet fluid leads to highly 
relativistic random motions within the medium and as a result the 
radiation from each sub-region is narrowly beamed in the direction 
of the local motion. This model is characteristic of a class of mod- 
els in which the observed variability time scale is determined by a 
process that takes place in the frame of the relativistically moving 
sub-region. Thus, both the observed variability time scale and the 
luminosity are affected by the motion of the sub-region (in con- 
trast to the subjets model where only the luminosity is modified). 
These two models are quite similar in spirit, but they differ in de- 
tails which become important in the case of extreme objects such 
as PKS 2155-304. 

In the case of the subjets model, using our canonical obser- 
vational constraints from PKS 2155-304, viz., flare duration T ~ 
20000 s, variability time St ~ 300 s, pulse duty cycle ^ ~ 0.4, 
reconnection inflow speed /3' ~ 0.1, we find that the model re- 
quires the Lorentz factor of the subjets, as measured in the mean jet 
frame, to be 7^ ^ 2500. Given that 7^ ~ y/a, where a is the mag- 
netization parameter of the jet fluid, this value seems unreasonably 
large. Perhaps of even greater concern, the model requires the to- 
tal number of independent reconnection sites, or sub-regions, to be 
ritot ~ 10^. Both requirements become less stringent if we assume 
extreme values for the time scales, e.g., T ~ 4000 s (i.e., setting 
T equal to the duration of the observations, assuming that the flare 
started exactly when the observations began) and 5t » 600 s (the 
maximum duration of each pulse, see lAharonian et aLluOOTT) . For 
this choice, r/(5t ^ 7,^ ~ 0.8, and we find 7^ « 170, ntot ~ 10^. 

The crux of the problem in the subjects model is the rela- 
tively low fiduci al value we assume for /3' (based on the work of 
lLvubarskvll2005n . The velocity /3'c is the speed with which mag- 
netized fluid moves in towards the reconnection point. Since this 
velocity determines the variability time, a low value of /3' implies 
a correspondingly small size I' of the reconnection cell. To com- 
pensate, 7j has to be very large so that the small amount of energy 
available within one reconnection cell is strongly beamed to give 
the luminosity observed in a single pulse of the light curve. 

It is possible that the arguments leading to /?' ~ 0.1 could be 
circumvented, allowing a reconnection inflow speed close to c. It is 
also possible that some (unknown) mechanism other than reconnec- 
tion produces the subjets and that this mechanism transports energy 
to the subjets at the speed of light. However, even if we set /3' — 1, 
this is still only marginally acceptable. Using fiducial values for the 
other parameters we find 7^ « 80, ntot ~ 10^. We obtain a rea- 
sonable solution only if we choose both a high value of /?' ~ 1 and 
extreme values of the time scales, T ~ 4000 s, St ~ 600 s. In this 
case, we find 7^ « 5, ntot ~ 10"^. Or, we could choose /?' ~ 1 and 
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instead of changing T and 5t, we could select a smaller BH mass, 
Mg « 0.5, which gives 7^ ^ 20, ntot ~ 4 x lO". 

We thus conclude that the subjets model works only if mag- 
netic reconnection in relativistic plasmas proceeds at the speed of 
light, much faster than currently thought. This alone may not be 
sufficient to explain PKS 2155-304. We may also need to push the 
observational estimates of time scales or BH mass to extreme val- 
ues. 

The situation in the case of the turbulence model is quite dif- 
ferent. In contrast to the subjets model, here, even with fiducial 
values of parameters, the implied conditions in the jet are quite rea- 
sonable. For instance, with T ■^ 20000 s, 5t ^ 300 s, .^ ~ o.4, 
we find j'j ^ 5 — 7, ntot ~ 10'^. Such conditions appear quite 
reasonable for a relativistically moving high-energy source. 

The weakness of the turbulence model is that it is short on 
details. We do not have a physical model of what causes the tur- 
bulence and what the limits of this process are. (This is in con- 
trast to the subjets model where we have a very specific process 
in mind - reconnection - which immediately gives a physical con- 
straint /?' ~ 0.1.) In fact, the word "turbulence" itself is merely a 
code for random motions. The actual dynamics may not be truly 
turbulent in the sense it is normally understood. Certainly, we can- 
not have hydrodynamic turbulence since motions with 7^ ~ 5 
are highly supersonic with respect to the maximum allowed sound 
speed, Cs,max ~ c/-\/3. In principle, MHD turbulence is compati- 
ble with the model since wave speeds can reach up to Lorentz fac- 
tors Fwavc ~ \/^, where a is the magnetization parameter There- 
fore, with a sufficiently strongly magnetized medium {a 3> 1), 
random motions with 7^ « 5 — 7 could be supported. 

The heating of particles and the production of radiation is also 
unexplained in the turbulence model. In the subjets model, recon- 
nection automatically produces relativistic beams of particles and 
it is not hard to imagine that these particles will radiate in the 
ambient magnetic field. (Note, however, that the dissipation effi- 
ciency in reconnection may be rather low unless the reconnecting 
fields coming in from the two sides are practically parallel, see 
iLvubarskvl 120051) . In the case of the turbulence model, the most 
obvious candidate for particle heating is a shock, but this is un- 
likely to work sinc e highly mag netized shocks are very inefficient 
at particle heating dKennel & C oroniti 1984; Naravan et al. 2011). 
Thus, it may be necessary to invoke heating through reconnection 
inside the turbulent blobs (which would immediately introduce in- 
efficiency through a /?'-like factor, just as in the subjets model). 
Perhaps the most promising idea is that the particles are heated di- 
rectly by waves in the turbulent magnetized medium, but the details 
remain to be worked out. 

Summarizing, we find the situation far from clear. The subjets 
model with magnetic reconnection provides a natural way to pro- 
duce the required "jets within a jet." How ever, given our c urrent 
understanding of relativistic reconnection l lLvubarskvll2005r) , viz., 
that the inflow speed towards the reconnection site can be no larger 
than a tenth the speed of light, this model has great difficulty ex- 

i ilaining the observations in PKS 2155-304 . The turbulence model 
Naravan & KumanI 20091 : iLazar et alj|2009l) apparently has no dif- 
ficulty fitting the data. However, all we have is a broad outline of 
this model, and there is no physical picture of how the "turbulence" 
in this model is produced or how this turbulence heats particles and 
produces the observed radiation. 

In §4 we analyze the pair production opacity for TeV photons 
in the "jets in a jet" model. We show that collisions of beams emit- 
ted from different sub-regions are common and hence the opacity 
due to these beam-beam collisions is more important than the opac- 



ity within a single emitting sub-region. Allowing for this effect we 
find that, regardless of which version of the "jets in a jet" model we 
consider, the minimum bulk Lorentz factor Ft of the jet is ~ 25, 
which is lower than the limit ~ 50 obtained by iBegelman et al.l 
( 120081) . Although the reduction in the value of Ft is only a factor 
of 2, it is nevertheless a significant revision since the new value 
is closer to the typical bulk Lorentz factors ~ 10 — 20 found in 
relativistic jets in other blazars. 
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APPENDIX A: BAYESIAN ESTIMATE OF THE 
DURATION OF THE FLARE IN PKS 2155-304 

Let us suppose that an astronomical source is on when we first start 
observing it. It remains on for a time Tobs and then shuts off. We 
would like to estimate the total on-time T of the source, including 
the time it was on before we began observing. 

Let us suppose that the source was observed previously a time 
Tmax before the current turn-off time and that, at that time, the 
source was off. Therefore, the minimum and maximum possible 
on-time of the source are Tobs and Tmax, respectively. We wish to 
evaluate the probability distribution of T between these two limits. 

In the absence of any other information, the probability distri- 
bution of T is logarithmically flat, i.e.. 



P{T)dT oc dlnT = dT/T. 



(A.1) 



For a given value of T, the chance that the observed on-time will lie 
between Tobs and Tobs + dTohs is clearly uniform for all values of 
Tobs between and T. Thus, for a true on-time T, the probability 
of an observed on-time Tobs is simply 



P(robs|r)drobs = dTobs/r, o < Tobs < r. 



(A.2) 



This distribution is normalized such that the total integrated proba- 
bility is unity. 

Bayes' theorem states that P(r|Tobs). the probability that the 
true duration is T, given an observed duration Tobs, is 

P(r|Tobs) = CP(robs|r)P(r), (A.3) 

where C is a normalization constant. Thus 

P(T|robB)dr = (robs/T')dr, Tobs < T < Tn^ax, (A.4) 

where the term Tobs in the numerator on the right-hand side takes 
care of the normalization (assuming for simplicity that Tmax 3> 
Tobs). The probability distribution P(T|robs) peaks at T = Tobs 
(which is reasonable), but it has a long tail extending all the way to 
T — Tmax. The median and mean values of T are easily calculated 
from this probability distribution: 



2To, 



(A.5) 



J mean — tn(_/ max/ -^ obs) -^ obs . V^-O) 

For the particular flare of interest in PKS 2155-304, we have 
Tobs ~ 4000 s, so Tmcdian ~ 8000 s. To compute the mean we 
need to know Tmax. If the source was observed the night previ- 
ous to the flare and it was off at that time, then Tmax ~ 10^ s, 
and we obtain Tnean ~ 3Tobs ~ 12000 s. However, if the pre- 
vious observation in th e off state was as long ago as Sept. 2003 
JAharonian erai]|2005h , then Tncan ~ lOTobs ~ 40000 s. 
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Table Bl. Comparison of notations 





This work 


Narayan & 
Kumar (2009) 


Lazar 
et al (2009) 


Giannios 
et al (2009) 


Radial distance from the BH 


R 


R 


R, Ro 


R = rRg 


Overall flare duration 


T 


t-burst 


T 


- 


Individual pulse time scale 


St 


Cvar 


St 


*/ 


Duty cycle 


€ 


Si 1 


rip 


- 


Bulk Lorentz factor 


Tf, 


r 


r 


r, 


(Flare duration)/(ii;/2r^c) 


fd 


1 


d 


- 


Random Lorentz factor 


^^ 


7t 


7 


Tco 


Reconnection region size 


I' 


- 


/' = V'SjR 


I 


Turbulent blob size 


V 


Te 


I' = iiR 


- 


Solid angle of each beam 


^'j^^fjhf 


^ht 


^^ 


i/4r?o 



APPENDIX B: COMPARISON OF NOTATIONS IN 
DIFFERENT WORKS 

Since the present paper dr a ws on a number of e arlier works, notably 
iNaravan & Kumaj O009h . iLazar et all ( |2009|) and iGiannios et all 
( l2009f) . we compare tlie notations in Table 1 . 



Sari R., Piran T, 1997, ApJ., 485, 270 

Tchekliovskoy A., McKinney J. C, Narayan R., 2008, MNRAS, 

388,551 
Volpe R, Rieger F. M., 201 1, ArXiv e-prints 
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